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Abstract 

We investigate the bound states of heavy mesons with finite masses to a 
classical soliton solution in the Skyrme model. For a given model Lagrangian 
we solve the equations of motion exactly so that the heavy vector mesons are 
treated on the same footing as the heavy pseudoscalar mesons. All the energy 
levels of higher grand spin states as well as the ground state are given over a 
wide range of the heavy meson masses. We also examine the validity of the 
approximations used in the literatures. The recoil effect of finite mass soliton 
is naively estimated. 
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I. INTRODUCTION 



The bound state approach of Callan and Klebanov has widened the apphcabihty of 
the Skyrme model up to heavy flavored baryons. In this approach, heavy baryons are 
described by bound states of the soliton of an SU (2) chiral Lagrangian and the heavy mesons 
containing the corresponding heavy flavor. This picture originally introduced for the strange 
baryons was suggested to be applied for the heavy baryons containing one or more heavy 
quarks such as charm (c) and bottom (b) quarks by Rho, Riska, and Scoccola 03 . The 
results on the mass spectra 0] and magnetic moments [§] for the charmed baryons were 
found to be strikingly close to the quark model description. 

Although qualitatively successful, when straightforwardly extended to the heavy flavors, 
the way of treating the heavy vector mesons in the traditional bound state approach |^ 
cannot be compatible with the heavy quark symmetry 0. As far as the strange flavor is 
concerned, on the analogy of the case of p and vr, the vector mesons K* may be integrated out 
via an ansatz in favor of a combination of the background and the pseudoscalar meson field 
K. (See Sec. Ill for its detailed form.) However, such an approximation is valid only when 
the vector meson mass is much larger than that of the pseudoscalar meson. Furthermore, 
the ansatz suppresses the vector meson fields by a factor inversely proportional to the vector 
meson mass, while the heavy quark symmetry implies that the role of the vector mesons 
becomes as important as that of the pseudoscalar mesons in the heavy quark mass limit. 

In the work of Jenkins, Manohar, and Wise ^ followed by a burst of publications P|-p^ 



this problem has been neatly solved out. There, the bound state approach is applied to the 
heavy meson effective Lagrangian that explicitly incorporates the heavy quark symmetry 
and the chiral symmetry on the same footing. The resulting baryon mass spectra show the 
correct hyperfine splittings consistent with the heavy quark symmetry. It emphasizes the 
essential role of the heavy vector mesons in the binding mechanism. This model is thoroughly 
studied in Ref. |]14| for higher spin states and is further applied to the pentaquark exotic 
baryons |]15 . 



However, these works are carried out in the limit of both the number of color Nc and 
the heavy quark mass mg going to infinity, where the soliton and the heavy mesons are 
infinitely heavy and so sit on the top of each other. This naive picture would definitely 
yield larger binding energies to the heavy-meson-soliton system. Besides, degeneracies and 



parity doubling in the heavy baryon spectroscopy of Ref. [|T^ would be just an artifact 



originating from this assumption. In Ref. |16], the kinetic effects of finite heavy mesons 



have been estimated up to 1/mQ order for the ground state. It is shown that the kinetic 
effects amounts to about 0.3 GeV in the charm sector, which is not small compared with the 
leading order binding energy of 0.8 GeV. Such a correction is expected to be more serious 
for the loosely bound pentaquark exotic states with a binding energy about 0.3 GeV. 

In this paper, on our way of investigating the pentaquark exotic baryons and the heavy 
baryon excited states in a more realistic way we generalize the work of Ref. [|16| to 



obtain the heavy meson (and antiflavored heavy meson) bound states of higher grand spin 
and radially excited bound states. For a given model Lagrangian, we solve the equations 
of motion exactly without adopting any approximations. It enables us to investigate the 
bound states over a wide range of the heavy meson masses from the strangeness sector to 
the bottom sector, and allows us to examine the validity of the approximations used in the 
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previous works. We work simply in the soUton- fixed frame, neglecting any recoil effects due 
to the finite soliton mass. At the end, a naive estimation of the soliton recoil effects is made 
by replacing the heavy meson masses by their reduced ones. 

This paper is organized as follows. In the next section we briefly describe our model 
Lagrangian for completeness. Then, in Sec. Ill, we derive the equations of motion for the 
heavy mesons under the influence of the static potentials and analyze the symmetries of 
the equations. In Sec. IV we discuss the spherical solutions of the free Proca equations 
before we give our numerical results for the wave functions of the bound states in Sec. V 
where the energy levels of the soliton-heavy-meson system are presented and discussed. 
We summarize the detailed formulas for the grand spin eigenstates in Appendix A and the 
equations of motion for the radial functions in Appendix B. 



II. MODEL LAGRANGIAN 

We work with a simple model Lagrangian for a system of Goldstone bosons and the 
heavy mesons, which possesses the chiral symmetry explicitly and restores the heavy quark 
symmetry as the heavy meson mass goes to infinity. 

As for the Goldstone bosons, we adopt the Skyrme model Lagrangian 

£SM ^ llt.^d^u^d^u) + ^tT[U^d,U, U^dM?, (2.1a) 

where /^r is the pion decay constant (~ 93 MeV empirically) and U is an SU (2) matrix of 
the chiral field; viz., 

U = e'^/^% (2.1b) 
with M being a 2 x 2 matrix of the pion triplet 

Here, the chiral SU{2)l x SU{2)ji symmetry is realized nonlinearly under the transformation 
of U: 

U — > LUR\ (2. Id) 



with L G SU{2)l and R G SU{2)ji. By the help of the Skyrme term with the Skyrme 

'SM 

■-M 



parameter e, the Lagrangian supports a stable winding- number- 1 soliton solution under 



the "hedgehog" configuration 

[/o(r) = exp[ir ■ rF(r)], (2.2) 

where the profile of F(r) is subject to the boundary condition, -F(O) = vr and F(r) — as 
r becomes infinity. The solution carries a finite mass of order Nc and the winding number 
is interpreted as the baryon number. The pion decay constant and the Skyrme parameter 



3 



are fixed as /7r=64.5 MeV and e=5.45, respectively, so that the quantized soliton fits the 



nucleon and A masses 18 



Now, consider j'^ = 0~ and 1~ heavy mesons with quantum numbers of a heavy quark 
Q and a light antiquark q. We take them as point-like objects described by the fields $ and 
$*, respectively, which form SU{2) anti-doublets. For example, if the heavy quark is charm 
flavored, they are the D-meson anti-doublets: 

$ = and $* = (2.3) 

Their conventional free field Lagrangian density is given by 

= a^$a'^$t _ ^|$<|,t _ 1$*M^<|.4 + ml,^*''^*^^ (2.4) 

where = 9^$* — d^^*^ is the field strength tensor of the heavy vector meson fields 
and m$ (m$,) is the mass of the heavy pseudoscalar (vector) mesons. 

In order to construct a chirally invariant Lagrangian containing $* , and their couplings 
to the Goldstone bosons, we need to assign a transformation rule with respect to the chiral 
symmetry group SU{2)l x SU{2)r to the heavy meson fields. For this end, we introduce 

e = t/^ (2.5) 

which transforms under the SU{2)l x SU{2)ji as 

e ^ r = L^^^ = ^^R\ (2.6) 



where t9 is a local unitary matrix depending on L, R, and the Goldstone fields M{x). Then, 
we can construct a vector field and an axial vector field as 

The vector field V"^ behaves as a gauge field under the local chiral transformation while the 
axial vector field transforms covariantly; viz.. 



(2.8) 



In terms of the chiral transformations of $ and $* fields are expressed as 

$ ^ $' = $i9t and ^ = <^l§\ (2.9) 
Therefore, a covariant derivative can be defined as 



D,<^1 ^ ^l(d, +v;t). 



(2.10) 



Such a complication as the multiplication of the field from the right hand side is due to 
the anti-doublet structure of the heavy meson fields $ and $* . 
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Given the above definitions, one can write down the chirally invariant Lagrangian for $ 
and $* with their couplings to the Goldstone bosons. Up to the first derivatives acting on 
the Goldstone boson fields, it has the form W9 



(2.11) 



where £0123 = +1, fq and gq are the $*$M and <I>*<I>*M coupling constants, respectively, 
and the field strength tensor is now newly defined in terms of the covariant derivative of Eq. 

( prmi ) as 

= - (2.12) 

In this work, instead of restricting the heavy meson masses, m$ and m$., to their 
experimental values, we will take them as free parameters given by a formula 



3k 

m$ = m$ — , 

K 



(2.13) 



4m<i, ■ 



with the averaged meson mass m$ = |(m$ + 3m$*). Given in Fig. 1 are the heavy meson 
masses evaluated by the formula ( p.l3|) with the parameter k fixed by the charmed meson 
masses rriD and rriD*. One can see that the formula reproduces the heavy meson masses 
remarkably well not only in the bottom sector but also in the strangeness sector. 

As for the coupling constants, since there is no experimental data available for them 
(except the upper limit Igql"^ < 0.5 estimated via the D* decay width and D*~^ decay 
branching ratio), we will adopt the heavy-quark-symmetric relation [|l^ 

fQ = 2m^,gQ, (2.14) 

and take the nonrelativistic quark model estimation for gq^ 

gQ = -0.75 (2.15) 

over a wide range of heavy meson masses. 

III. EQUATIONS OF MOTION FOR HEAVY MESON EIGENMODES 

What interests us is the heavy meson eigenstates bound to the static potentials formed 
by the fields, 

V^ = {V',y) = iO,-tiTXv)vir)), 

A'' = {A^, A) = (0, i[ai(r)r + a2{r)rT ■ r]), ^""''^ 
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with v{r) = [sm^(-F/2)]/r, ai(r) = (sinF)/r, and a2{r) = F' — (sinF)/r, (Hereafter, we will 
denote the derivative with respect to r and t by a prime and a dot, respectively. That is, 
/' = df/dr and / = d^f.) which are provided by the soliton configuration ( p.2| ). We shall 
work in the soliton-fixed frame, neglecting any recoil effects of finite mass soliton. In Sec. 
V, we will roughly estimate such effects. 

The equations of motion for the heavy mesons can be read off from the Lagrangian (|2.11|) 

as 

D^D>^^^ + m|$t = fQA^^;\ (3.2) 
D,,<l>*^^t + ^|,$-t = -fQA''<^^ + gQe^^'^^PAx<l>;l. (3.3) 

In order to avoid any unnecessary complications originating from the anti-doublet structure 
of $ and $*, we will work with and instead. 

An auxiliary condition corresponding to the Lorentz condition 5^$* = for the free 
Proca equations can be obtained by taking a covariant divergence of Eq. (|3.3|): 



m|,Z},<l>-t = _Z)^Z)^$*M-t _ J-qDM"'^^) + gQe^'^'M^x'^t)- (3-4) 

Note that none of the terms in the right hand side vanishes identically and that the auxiliary 
condition, thus, cannot be simply reduced to D^^*^ = 0. Thus, different from the free Proca 
equations for the spin-1 field, it is very difficult to eliminate, for example, $q by using Eq. 
( p.4|) at this level. Actually, such an elimination is not an indispensable process in solving 
out the Proca equations as we shall see in the next section. 

In terms of = ($*'^, $*) and using the form of (|3.1| ), Eqs. ( p.2[ - p.4|) can be rewritten 
explicitly as 

$t _ D2^t + ^2 = _ j^A ■ (3.5) 
(D^ - m|.)$S^ = D ■ $*t _ 2gQA ■ (D x $*t)^ (3.6) 

+ D X (D X $*t) + m|.$*"f = B^l^ - fQA<^^ - 2^qA x ($*+ - D$o^), (3.7) 
m|,(<i>;^ - D • $*t) - D ■ [D X (D X 

= /qD ■ (A<l>t) + 2(7q[D ■ A X ($*^ - D$;t) + A ■ D X ^% 



(3.8) 



with D = -V + V. 

To solve them, we need to know the symmetries of the equations of motion. First of 
all, they are invariant under parity operations under which the heavy meson fields and the 
static soliton field transform as 

^[r,t) — > — $(— r,t), 
mr,t) ^ +n{-r,t), $*(r,t) ^ -$*(-r,t), (3.9) 
f/(r, t) — . t/t(-r, t), (V, A — . -V, +A), 

where we have used that the heavy mesons have negative intrinsic parity. Next, due to the 



correlation of the isospin and angular momentum in the hedgehog configuration ( |2.2| ) and 
consequently in the static fields ( |3.1|) , the equations are only invariant under the "grand 
spin" rotation generated by the operator 
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K = J + I = L + S + I, 



(3.10) 



where L, S, and I, respectively, denote the orbital angular momentum, spin, and isospin op- 
erator of the heavy mesons. (See Appendix A for their explicit forms and the corresponding 
eigenbases.) Thus, the eigenstates are classified by the grand spin quantum numbers {k, k^) 
and the parity vr. 

For a given grand spin (fc, k^) with parity vr = (— l)'^^^/^, the general wavefunction of an 
energy eigenmode can be written as 

$t(r,t) = e+-V(r)3^g(f-), 

$;t(r,t) = e+^-VSW3^£(i-), (3-11) 
**t(r,t) = e+-* [pl{r) r3^(g(f) + ^l{r) L3^g(r) + ^l{r) G3^g(r)" 

with five radial functions ({>{r-) and v'a(r) (a = 0, 1, 2, 3). Here, G is an operator defined 
as G = — i(r x L), and y^l{i) are the spinor spherical harmonics obtained by combining 
the orbital angular momentum eigenstates and the isospin eigenstates. (See Appendix A for 
further details.) Note the anomalous sign convention of the energy in the exponent for the 
time evolution of the eigenmodes. It is due to the fact that we are working with and 
instead of $ and $*; thus, the eigenenergy oj is positive for the bound states of heavy mesons 
and is negative for the "antiflavored" heavy mesons of Qq structure. Substituting Eq. (p.ll|) 



into Eqs. ( |3.5| - |3.7D , we can obtain the equations of motion for the radial wavefunctions. Their 
explicit forms are given in Appendix B. 

In the literatures, in order to reduce the complexity of the equations of motion, the 
solutions have been found approximately by adopting proper ans"' atze. In case of sufficiently 
heavy mesons, one may drop the l/mq and higher order terms in Eq. ( p.4|) which leads us 



to a simple auxiliary condition |]T6 



=0 (Ansatzl). (3.12) 

It enables one to eliminate $q easily in favor of the other three fields In order to be 
consistent, one should also drop the higher order terms in the equations of motion for v^i 2 3- 
On the other hand, in the limit of light mesons one may use the fact that the vector 
mesons are much heavier than the pseudoscalar mesons. Then, the most dominant terms 
in the Lagrangian are those for the vector meson mass and the $$*M interactions with the 
coupling constant fq, which reads 

m|.<|.*^<|.*t + /Q(<|.A^$*t + A^$t) 



with e abbreviating 2gQ/mq,*. Here, we have used the relation (|2.14| ) for the coupling 



constant fq. Equation ( p.l3| ) suggests us another ansatz for the vector mesons 0: 

$*t = _^^^$t (Ansatz II). (3.14) 
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In terms of the radial functions, it reads 

if*o{r) = 0, ^l{r) = le{ai + a2Mr), ^;{r) = ^e^^aiif{r) = -</?^(r). (3.15) 

However, the equation of motion for the field $ (or (f) should not be the one obtained simply 
by substituting Eq. (|3.14| ) or ( p.l5| ) into the corresponding equation to eliminate $* . Under 



this ansatz, the Lagrangian ( p.ll| ) is reduced to that of the pseudoscalar field $ only: 

(3.16) 



-Agl^A^A^^^^ + gQe''e^^^P{D^<^A,Ay^A,<^^ + <^ApAxA,D^<^^), 



which yields the equation of motion for $ as 

(1 + e'^A ■ A)<l>t + 2^q£2a • (A x A)$t _ D^^t + (m| - Ag^A ■ A)<l>t 
-e2(A X D) • (A X D)<l>t = 0. 



(3.17) 



Substituting $t(r,t) = e*'^V(^)3^^(r), we can obtain the equation of motion for ^{r) as 

h{r)^" + h'^' + [f{r)u^ + X{r)uj - m| - V,s{r)]^ = 0, (3.18) 

where 

h{r) = 1 + le^al h\r) = ^[r^h{r)] , 
f{r) = l + \e%a, + a2f + 2af^, 
A(r) = \£^gQal{ai + 02), 



<j ( ^ + 2t;j af + ai(a; + 4) + 2t;j + (ai + 02)' I fa? + 



(3.19) 



The similarity of the equation to those of Refs. is remarkable. Those e-dependent terms 
play an important role in obtaining the bound states in the light meson mass limit. Note 



that the factor e is not necessarily small for the Ansatz II to be a good approximation [|2T |. 
In the strangeness sector, for example, we have e ~ —0.85 [in (e/^)~^ unit]. 

However, as the mass of the involving heavy flavor increases, the mass of the pseudoscalar 
mesons becomes comparable to that of the vector mesons, and Ansatz II cannot be expected 
to work well. Note that it suppresses the role of the vector mesons by the factor e, while 
the heavy quark symmetry implies that they are not distinguishable from the pseudoscalars 
as far as the low energy strong interactions are concerned. In Sec. V, we will give the exact 
solutions with the approximated ones using ( |3.12| ) and ( |3.14D . 



8 



IV. SOLUTIONS TO FREE EQUATIONS 



Before proceeding, we digress for a while to get an insight from the spherical wave 
solutions to the free equations, which are obtained by inserting (|3.11|) into C^^^'^ of 



A 



if" + -if' - (m| - Lu'^ + '-^)(p = 0, 



A. 



2 .2 , * , „*/ I "^T I *i ^ 



= ^Vl + — (^3 + -<^3)) 



9 1 1 

*//, */ / 2 2\ * */ * 



and the Lorentz condition 



2 A 

^<^0 = + -Vl) - —^3- 



(4.1) 
(4.2) 
(4.3) 
(4.4) 
(4.5) 

(4.6) 



For a comparison, we are insisting to find the solutions in the form of Eq. (|3.11|) . In case of 
free fields, ip{r) and (flir) are completely decoupled to the other radial functions and Eqs. 
(4.1) and (|4.4| ) lead to the solutions as 



V{r)=j.{qr) and ^l{r) = j.{q*r), 



(4.7) 



where jg^{x) is the spherical Bessel function of order i±{= k ^1/2) with = — m| 
and g*^ = ~ m|. . That is, we have found two of four independent solution sets of the 
equations: 



Solution I : 



Solution II : 



<^>nr,t) = e^"*j,Jgr)#j^(r), 

f $t(r,t)=0, 
*S^(r,t) = 0, 

$*t(r,t) = e*-*j (g*r) Ly, 



k nih 



(4.8) 

(4.9) 



up to some proper normalizations. Note that there is no conjugate field of $q in the "Solution 
II," which is consistent with the Lorentz condition ([4.6|) . 

The other two solution sets can be found by the following way. With the help of the 
Lorentz condition (|4.6|), Eqs. ([4.2|), ( |4.4| ), and ( |4.5| ) can be rewritten as 



r 



r 
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*2 



(4.10) 
(4.11) 
(4.12) 



By diagonalizing the v^i-y^s mixing part as 



we can obtain the two remaining 

r Voir 

Solution III <^l{r 

i Voir 
Solution IV : <^ ^1 

for the TT = (—1)^+1/2 states and 

r Vl{r 
Solution III A ^\ 
\ Vl 




4(4 + i: 



Vl 



independent solutions as 

= ik + l/2){q*/u;)j,^,/,iq*r), 
= {k+l/2)j,^,^,iq*r), 
= 3k-i/2{(fr), 

= {k + 3/2){q*/u)j,^,/,{q*r), 
= (A; + 3/2)j,+3/2(gV), 
= -Jk+3/2i(l*r), 



r 



Solution IV : ■ 



,,., = ik + l/2){q*/u;)j,_y,{q*r), 
{r) = {k + l/2)j,^,j,{q*r), 
= -Jk+i/2i(l*r), 

,,., = {k-l/2){q*/uj)j,_,^,{q*r), 
(r) = (/.-l/2)j,_3/2(gV), 
'Plir)=Jk^3/2i(l*r), 



V- 



for the 71 = states. 

From these free solutions we can learn the followings: 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



1. In principle, we can obtain the solutions in the form of Eq. ( |3.11| ) without eliminating 



one of $* at the level of Eqs. ( p.5H3.8| ). It is enough to keep in mind that, due to the 
dependence in field variables (and consequently in the equations of motion), we can 
obtain only four independent solutions. 

2. The ^l-^l mixing is due to the fact that our grand spin eigenstates, ^3^1^^ ^^3^^^^, 
are not the eigenstates of the orbital angular momentum L^. If we have used iy^mj,, 
the conjugate radial functions would have been decoupled. For example, the solution 
III of Eq. (f4.14|) can be written in a form of 



<I>t(r,t) = 0, 



Solution III: <! %\r,t)=iik + l/2)iq*/u)j,^,Jq*r)yi^J^^iv 



#*t(r,t) 



Jk-i/2iq*r) [(A; + l/2)r3^^(r) +G3;^^^> 



(4.18) 



{2k + l){k+l)j,_,/,{q*r)y 



( 



3. The special form of Eq. ( |4.3| ) and ( [B3| ) with no derivative on Lfl is also ascribed to 
our choice of the grand spin eigenstates. Although one can easily express ip^ in terms 
of the other fields by using the equation 
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= c^-mi-A^/H ' ('-'^^ 

it cannot be used in eliminating Lp\ from the equations. Let's see what happens if 
one proceeds along this direction. Substituting Eq. (|4.19| ) and its derivative into Eqs. 
121) and fO) leads to 



m: 



cij^ — m|, — Azp/r^ ° ti;2 — m|, — Azp/r^ ^ ' 



(4.20) 



where po,3 are the terms without containing the second order derivatives. In order to 
obtain (/jQ'g, we should solve the above linear equations. However, the matrix of the 
coefficients is singular at the origin and, thus, one cannot obtain the correct solutions 
in this way. 



V. NUMERICAL SOLUTIONS AND DISCUSSIONS 



In solving equations numerically, there is a considerable freedom. We will solve Eqs. ( plf) 



and (P^B6| ) by taking yj', ^Jq, yj^, yjj, V's' ^"^^ eight independent variables. Note 

that Eqs. (pi]), (pif), and ( p5D are the second order differential equations for yj, yjj, and ^93, 
while Eqs. ( |B3[ ) and ( |B6D are the first order ones for ip^ and y?^, respectively. On the other 
hand, we can take Eqs. ( piP , (p2D, (pi]), and (|B5|) as the second order differential equations 
for (y9Q, (^25 and and use Eqs. ( P^D and ( PB] ) to eliminate (p\ and y^f appearing in those 
equations. One may be tempted to eliminate ip\ by using only Eq. ( |B3| ), since it does not 
contain any derivative on ip\. However, the numerical program according to such a recipe 
meets serious instabilities near the origin due to the singularity in a linear equation similar 
toEq. i WM )- 

In Fig. 2, we give the radial wavefunctions for a few eigenstates (solid lines). For a 
comparison, we present the approximate solutions obtained by using Ansatz I (dashed lines) 
and Ansatz II (dash-dotted lines) as well. In case of the Ansatz II, ip\ 2 are obtained from 
the relations (|3.15| ) with the solution of Eq. ( p.l8| ) substituted for ip. As we have expected, 
Ansatz II works only in the limit case; the approximate radial functions of Ansatz II are 
quite close to the exact ones in the strangeness sector but it becomes a bad approximation in 
the charm sector. One can see that the factor e of the Ansatz II suppresses the vector fields 
too much. On the other hand, Ansatz I works quite well both in the strangeness sector and 
in the charm sector. In the bottom sector, the approximate solutions are indistinguishable 
from the exact ones. This observation supports the approximation of Ref. ||16|| . 



From Fig. 3, one may arrive at the same conclusion for the Ans"'atze I and II, where we 
present the eigenenergies of the k'^=\'^ ground state, \^ radially excited state, and \ state 
as functions of the heavy pseudoscalar meson mass m<i,. That is, Ansatz I works well over 
a rather wide range of the heavy meson masses while Ansatz II works well only in the limit 
cases, say, m$ < 0.8 GeV. Therefore, the Ansatz II can be justified only in the strangeness 
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sector as in Ref. Note that the Ansatz I and II could yield lower eigenenergies than 
the exact ones. However, it does not contradict our common sense that the exact solutions 
should have the lower eigenenergy than the approximate ones. Through the ans" atze we have 
altered the interactions of the heavy mesons with the soliton more or less. Note also that 
the equations admit one bound state in the strangeness sector without the higher derivative 
terms. In case of Ansatz II, the e terms are essential for the existence of the bound state. 

Given in Fig. 4 are the binding energies uj — of the bound states as functions of the 
heavy pseudoscalar meson mass m$. As the heavy meson mass increases, the fc'^ = | ground 
state energy decreases and it may reach down to the infinite mass limit —^gqF'ifd) ~ 0.8 
GeV. However, the kinetic effect turns out to be quite large and the binding energy does not 
come to the infinite mass limit value even when the heavy meson mass is increased up to 10 
GeV. Note that, as we have expected from the heavy quark spin symmetry, the eigenstates 
of A; = fcf ±1/2 become degenerate as the heavy meson mass increases. A special care should 



be taken in the quantization of these degenerate bound states. WM In Fig. 4, there appear 
many bound states for the heavy mesons in the charm and bottom sector. However, remind 
that we have neglected any recoil effects due to the finite soliton mass. At this point, we 
may naively estimate the effects of the soliton motion by replacing the heavy meson masses 
with their corresponding "reduced" masses defined by 

= ■ , (5.1 

m$ + rrtsoi 

where mgoi is the soliton mass. With mgoi = 867 MeV, we have /i^) ~ 590 MeV and fiB ~ 743 
MeV. From Fig. 4, one can see that only one (or at most two) bound state (s) could survive 
when the recoil effects are incorporated. This correction deserves to be studied further. 
As discussed in Refs. |I^ , P3| , the equations admit the bound state solutions with negative 



eigenenergies, which can be interpreted as the bound states of antifiavored heavy mesons. 
In Fig. 5, we present the radial functions of a few low-lying bound states in case of Q = b. 
It is apparent that the radial functions spread over a wider range of r than those of Q = b 
do. The binding energies of antifiavored heavy mesons are given in Fig. 6 as functions of the 
heavy meson mass. In this case, the energy change due to the recoil effects is comparable to 
the binding energy and the recoil effects seem to be crucial for the existence of the bound 
state. However, a naive estimation using ( p.l| ) shows that there still survive a bound state 
near the threshold, which leaves a possibility of stable pentaquark baryons. 

As a summary, we obtained the energy levels of the soliton-heavy-meson bound states. 
The equations of motion are solved exactly in a given model Lagrangian for the excited 
states as well as for the ground state. The calculation was also made for the pentaquark 
exotic baryons. However, to obtain the real mass spectrum of heavy baryons, we should go 
one step further; i.e., the system should be quantized for describing baryons of definite spin 
and isospin quantum numbers. Work in this direction is in progress and will be reported 
elsewhere. 

Note added. After completion of this paper, we were aware of recent work of Schechter 
and Subbaraman |]23|. In estimating the kinetic effects of heavy mesons, the authors adopt 



an interesting approximation to the equations of motion. Corrections due to finite soliton 
mass are roughly calculated in the same way as in this paper. 
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APPENDIX A: GRAND SPIN EIGENSTATES 

In this Appendix, we construct the grand spin eigenstates, i.e., the angular part of the 
wavefunction by combining the orbital angular momentum (L), spin (S), and isospin 

(I) of the heavy mesons. We first combine the orbital angular momentum and spin to get 
the angular momentum (J) eigenstates, and then combine the isospin 

To do this, we first find the spin operator Si {i = 1,2,3) for the vector mesons and the 
corresponding eigenstates. The Lagrangian ( |2.4] ) is invariant under an infinitesimal Lorentz 
transformation 



^'M _ ^At _|_ 



<l>(x) <l>'(x') = <l>(x), (Al) 

where e'"' = -e"'' and {S'"')ai3 = g'^aO'^fs - g^^g^a with g^^ = diag(l, -1, -1, -1). It defines 
conserved angular momentum 

^ J (A2) 

where V^^ is the canonical energy-momentum tensor and n*"(= dC%^'^ /d^*^) is the momen- 
tum conjugate to the field . Here, the indices run from 1 to 3. The first part corresponds 
to the orbital angular momentum L and the second to the spin angular momentum S of the 
vector mesons. The latter can be rewritten neatly as 

d^r (n* X + $* X n*'^), (A3) 

which defines the corresponding quantum mechanical spin operator acting on the vector 
meson wavefunctions as 

S'i=ieiX, 5'2 = 2e2X, 5'3 = ie3X, (A4) 

where ei, 62, and 63 are unit vectors along the x, and z axis, respectively. Strictly 
speaking, although they satisfy the correct commutation relations [Si, Sj] = iSijkSk, and the 
square is an invariant of the group, it is not the relativistic spin operator that commutes with 
all the generators of the Lorentz group. It is a good spin operator only in the rest frame, 
which is, however, enough for us to proceed. Actually, the invariance of the equations 
( P75| - ^.8[ ) under the grand spin rotation can be achieved only with this spin operator. 
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Eigenvectors of and S3 are found by taking suitable linear combinations of the unit 
vectors ei, 62, and 63 [^. Let's define 



e±i = T^(ei±2e2), ^^^^ 
eo = h, 

which satisfy 

§2 g = 2 e 

" with g = ±1,0. (A6) 

Then, the eigenvectors of and J3 for the pseudoscalar mesons are simply the spherical 
harmonics Yj>„^. (f), and for the vector mesons the "vector spherical harmonics" 

Yjirn.ir) =J2Yem{r)eg{£mlq\jmj), (A7) 

m,q 

where {imlq\jmj) is the Clebsch-Gordan coefficient of adding the orbital angular momentum 
and the spin. As a result of the angular momentum addition rule, there can be three different 
kinds of vector spherical harmonics with a given (j, rrij ) : 



Yiim,(r) with parity tt = -(-1)^' , ^^^^ 
Yjj±im^(r) with parity TT = +(-1) J 



where the parity vr incorporates the intrinsic parity of the vector mesons. 

The vector spherical harmonics can be generated from the spherical harmonics Yimii") 
by making use of certain operators: ||2^ 

h{j + 1) 



where G = — i(f x L). It enables us to carry out somewhat tedious calculations involving 
vector spherical harmonics through elementary vector algebra without referring Clebsch- 
Gordan coefficients. 

To complete the job, we combine the isospin to these angular momentum eigenstates. The 
quantum mechanical isospin operator acting on the isodoublet (not anti-doublet) structure 
$^ and is given by the Pauli matrices as 

I = |t. (AlO) 

We denote the eigenstates of and I3 as 
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X+1 



0^ 



(All) 



which have the properties I^x±i = fXii hx±^ = ^\x±i- 

Now, for the pseudoscalar mesons, the grand spin eigenstates are obtained simply by the 
spinor spherical harmonics: 



(i) for k = e + l/2 and tt 

= + 



-1) 



k+l/2 



;(+) 



^ k ' 

fii) for k = 



2k 



' k — nik 



2k 



1/2 and tt = (-1) 



fe-1/2 



(A12) 



^-^ ( 



Ik - nik + l 
2{k + l) 



Ik + nik + l 
2{k + l) 



which are related with each other by 



-[r-i) y^t' 



(A13) 

As for the grand spin eigenstates of the vector mesons, coupling of the isospin eigenstates 
to the vector spherical harmonics (|A7|) leads to six different grand spin eigenstates 3^1* for 
a given (fc, ^3): 



y^kra^J) = ^j^i^m,ir) Xg Ui rrij ^q\k nik), 

nijq 



(A14) 



with k = ji± 1/2 and ji = ii, ii ± 1. We list possible ji and ii in Table |. By the help of 

L i] 

1 



Eq. ( |A9|) , we can rewrite them in terms of f}^!^^,' ^3^^^^, and G3^|^^, for example. 



y^k^ 

^ k rni^ 



<^ k mj. 



y'l. 



T 

^^kmi^i 



2^' + 1 + 1) 



2j + l 
2j + l 



(A15) 



with j = k + l/2. For a practical convenience, we take ry^^^, L3^^^^, and G3^|^^ as six 

independent grand spin eigenstates instead of ^V^tife- 

The general solution to Eqs. (|3.5|-|3.8|) for a given (A;, nik) can now be written as 



$"^(r,t) 



+iujt 



(-) 



(+) 



?(-) 



^k nii^ ^ 



*{+) 



(+) 



(A16) 
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We may decompose them into two solution sets of definite parity as in Eq. (|3.11|) . 



APPENDIX B: EQUATIONS OF MOTION FOR RADIAL FUNCTIONS 



In this Appendix, we give the exphcit forms of the equations of motion for the radial 
functions ip{r) and (/^^(r) (a = 0,1,2,3). Substitution of Eq. ( p.ll|) into the equations of 
motion (|3.5|-|3.8|) leads us to the following coupled differential equations. 



2 

II \ I I 

ip + -Lp + 

r 



2 2 / A± 2t; 2 

u - m$ - — H /i± + 2v 

\ r 



r 



r 



2 2t; 



r 



iff! I I 

¥^2 + -¥'2 + 
r 



— H /i± + 7±/i±t; 



^2 



V 



+ hfQail±'f + gq {l±aigi + (ai + 02) [7±fi'2 + (1 + 7±)fi'3]} 



(Bl) 



(B2) 



(B3) 



(B4) 



*// I ^ */ I 2 2 I 2 * 

i - 7±t^) '^X - 7± (^t^' + - /i±t; (^^ + 7±t;^ - u; (^^ - 7±t; ) v^o (B5) 

-|/Qai7±V5 + gq {-7±aifi'i + (ai + ^2) [(1 - 7±)fi'2 - l±g^\} , 

+21^ (^-^^)/l-(^^' + ^)[/^T/2 + /^±/3] 



' 2 



[ai + a2)(v5' + ^V^) + (fli + a2)V + + 2t;)ai(y9 



(B6) 



+2^0 



where /j and gi (i=l,2,3) are functions defined by 
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D X ^ i [y^ + l^yi^2, + h 



(±) 



k rrik 

and A±, and 7± are the numbers defined through 



(B7) 

(B8) 



(B9) 



which also give r • G3^;^^^ = A*±3^imfc- Exphcitly, /j and gi {i— 1,2,3) are given in terms of 
cpl (a=0,l,2,3) as 



= - 

/2 = + 
/3= + 



'A. 



r 



92 = UJV*2 + l±v^*o, 



^3 = a;(^3 + ( - - 7±z; ) (Pq, 



and A±, and 7± are written in terms of k as 



A+ = (A;-l/2)(A; + l/2), 
/i+ = A; - 1/2, 
7+ = /i+/A+ = l/(/c + l/2), 



A_ = (A; + l/2)(A; + 3/2), 



-(A; + 3/2), 



7- = /^-/A_ = 

One can further derive some useful relations between them: 



-1/(^ + 1/2). 



(BIO) 



(Bll) 



(B12) 



A± = A^ + 2/i:p + 2, /i± = -(/x^ + 2), 
7± = -7t' /^±7± + /^t7t = 2, A± + /x± = A^ + i^^. 
Near the origin, the equations of motion behave asymptotically as 



A 



*// I 2 ^/ A-j- j , ^ 2 , 



+^ I ~(2 + /Xt)<^I - l^±(p*2 - ^^^^^(/?; - (^(fl' + ^(fl 



A 



r r r 



(B13) 



(B14) 
(B15) 
(B16) 
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2 (A±+/i±)(l+7±) . 

V2 + -V2 -Q. 



— Vl + 

r 



7± 



= ^2 - ^ + ^7±</'o , 

* / */ 2 ^ /i± ^ (A=p+/x=p) ^1 
^V'o - + -V'l) + — ¥^2 + V's f 



3?7(52V5o - ^2 



*/ I (At^2) ^ 

/^±V52 H V2 

r 



(B17) 



(B18) 



(B19) 



where we have used that, near the origin, the function F{r) behaves as F{r) = vr + F'{0)r + 
|F"'(0)r^ + ■ ■ which imphes 



v(r) - + 0(r), ai(r) ~ -F'(0) +0(r2), asfr) ~ 2F'(0) + 0(r^ 
r 

In the expressions ( |B14| - |BT9| ), the constants t], 61, and 62 are defined by 



r7 = (?QF'(0), 5i 



2mL 



0(l/mQ), 52 



2mL 



0(l/m^). 



(B20) 



(B21) 



With F'{0) ~ 2e/7r and exphcit heavy meson masses, one can evaluate the constants as 
61 ~ 0.65, 0.29, 0.11 and 62 ~ 0.31, 0.06, 0.01 in case of Q = s, c, b, respectively. Note 
that, due to the vector potential V[~ ^(r x T)/r, near the origin] in the covariant derivative, 
the singular structure of the equations ( |B14| - PT9| ) is quite different from that of the free 
equations (|]l]-|]6|). Since these equations are not all independent, they yield only four 
independent asymptotic solutions that are finite at the origin. We list them in Table 0, 
where c* and {i = I, II, III, a = 0, 1, 2, 3) for vr = (— l)'=+i/2 states are constants satisfying 

{k + 1/2)4 -ik + 3/2)4 - (A; + 3/2)^^ - (u - r])d^ = 0, 

(uj - 3r]52)ci, - [{k + 5/2) + 52ik - l/2)]c\ + (k - 1/2) [1 - 252(A; + 3/2)]4 (B22) 
+ {k + 3/2) p - 1/2) + 62{k + 5/2)]4 + 26i{k + l)c^ = 0. 

Since we have only two equations for five unknowns, we need to fix three of them. For 



example, we take (4 2 3 — 0'"^^ — 1)' (^^^ — ^3 



0,4' 



1), and (c 



Jii 



0,4" 



1), 



which yield three solution sets given in Table |T|. 

On the other hand, at large r where heavy mesons become free from the interactions 
with the soliton, the equations of motion (pll - p6|) approach asymptotically to Eqs. ( |4.1[ - 
[4.6|) . Since we are interested in the bound states with u;^ < m|, m|*, instead of the spherical 
Bessel function je{x) in the free solutions, the asymptotic solutions are written in terms of 
the modified spherical Hankel functions k£{x) satisfying 



--11 2 ~, 
k^ H kf^ — 

X 



1 + 



x^ 



k^ = 0, 



(B23) 
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with recurrence relations 

h-ii^) - h+iix) = [{2i+l)/x]k^{x) 



(B24) 

e-iix) + {£ + l)h+iix) = (2£ + 



Listed in Table ^ are four independent asymptotic solutions which decay exponentially at 

4 - and Q* = Jm^, 



infinity, where Q = \ m\ — uj'^ and Q* = \ m\, — uj'^. 
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FIGURES 



FIG. 1. Heavy meson masses and the mass formula (|2.13| ). The pseudoscalar meson mass m$ 
and the vector meson mass m^* obtained using Eq. ( 2.13| ) are given by sohd and dashed hnes, 
respectively. 

FIG. 2. Radial wavefunctions of the heavy meson bound states in arbitrary scale. Solid lines 
are those obtained by solving the equations exactly, while dashed lines and dash-dotted lines are 
the approximate solutions of Ans' atze I and II, respectively. 

FIG. 3. Binding energies, lo — m$, of k'^=^^ ground state, excited state, and ^ state. 
Solid lines represent the exact solutions and dashed (dash-dotted) lines correspond to Ansatz I 
(II). 

FIG. 4. Binding energies, u — of the bound states with k'^ as functions of the heavy meson 
mass. Solid (dashed) lines denote the positive (negative) parity states. 



FIG. 5. Radial wavefunctions of the antiflavored heavy meson bound states in arbitrary scale. 

FIG. 6. Binding energies, — m<j>, of the bound states as functions of the antiflavored heavy 
meson mass. Solid (dashed) lines denote the positive (negative) parity states. 
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TABLES 



TABLE L Eight grand spin eigenstates. 





TT = 


(_l)fc+l/2 






TT = 


(_l)fc-l/2 




state 


S 


e 


3 


state 


S 


i 


i 







3 


A; - 1/2 







3 


A; + 1/2 




1 


3 


A; - 1/2 




1 


3 


A; + 1/2 




1 


i-1 


A;+ 1/2 




1 


i-1 


A; - 1/2 




1 


i + 1 


A;+ 1/2 


^k ruk 


1 


i + 1 


A; - 1/2 



TABLE II. Four independent asymptotic solutions near the origin. 



parity 


TT = (-l)fc+V2 


TT = (-1)^-1/2 


Solution I 


^* ^ 4r'=-V2 + c>(r'=+3/2) 


^ ^ ^fc-1/2 + 0(rfe+3/2-) 

~ 0(r'=+V2) 

¥'1,2,3 ~ 0(r*^+3/2) 


Solution II 


(/?,^| ~ 0(r'=+5/2) 

iff, ~ 4V=-V2 ^ o(r'=+3/2) 

^* ^ r'=+V2 + OCr'^+s/^) 


~0(r*^+3/2) 
~ 0(r'=+V2) 

^ +^fc-l/2 _^ Q^^fc+3/2) 
^ _rfe-l/2 + 0(rfe+3/2) 


Solution III 


(p,^*2 ~ 0(r^+5/2) 

(/^S ~ 41^^-1/2 + 0(^fc+3/2) 
~ cI/V'=+V2 + (9(^fc+5/2j 
^* ^ j.k+1/2 ^ 0(^fe+5/2) 


<fl ~ 0(r^+3/2) 

~ 0(rfc+V2) 
(/?^ ~ (A; - l/2)r*^-V2 + o(r'=+3/2) 
^* ^ + o(r'=+3/2) 


Solution IV 


~ 0(r*^+V2) 
(p* ^ r*^-i/2 + 0(r'=+3/2) 
(/p^ ~ A+r*^-3/2 ^ (9(^fc+i/2) 

~ r*=-3/2 + (9(^fc+l/2^ 

~ //+r*=-3/2 + 0(r^+i/2) 


^ ~ 0(r'=+3/2) 

~ 0(r'=+3/2) 
(p* ~ 0(r'=+3/2) 
~ 0(r'=+3/2) 
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TABLE III. Four independent asymptotic solutions at large r. 



parity 


= (_l)fe+V2 


vr = (-1)'=-V2 


Solution I 


<1,2,3 ~ 


'P^h+i/2{Qr) 

V'(},1,2,3 ~ 


Solution II 


f*2 ~ h-i/2{Q*r) 
</^0,l,3 ~ 


^*2 ~ h+i/2iQ*r) 


Solution III 


(^* ~ 

^l^{k + l/2){Q*/u)h+i/2{Q*r) 

\,K + -L/^jfcfe-l/2W ^) 

</33 ~ ~kk-il2{Q*r) 


~ + l/2)(Q7a;)^,_i/2(Q*r) 
993 h+i/2{Q*r) 


Solution IV 


(^,(^* - 

^l^{k + ^/2){Q*/u)h+i/2{Q*r) 
(^t ~(fc + 3/2)fcfc+3/2(Q*r) 
h+3/2kk+3/2{Q*r) 


y?^ ~ 

y,* ^ - l/2)(Q7a;)^,_i/2(Q*r) 

~(fc-l/2)^,,_3/2 (Q*r) 
</'3 ~ h-z/2{Q*r) 
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